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coefficients  of  the  composite  polynomial  for  ice  velocity  at  a  point 
parameter  groupings  in  the  equation  of  motion  for  the  ice 
acceleration  of  the  ice 
river  width 

coefficients  in  the  recurrence  relation  forothogonal  polynomials 
speed  of  the  breaking  front 
dynamic  wave  speed  in  the  downstream  direction 
dimensionless  Chezy  conveyance  coefficient  of  the  ice 
apparent  or  added  mass  coefficient 

coefficients  in  a  linear  combination  of  orthogonal  polynomials  that  approximate 
the  ice  velocity  data 

ith  coefficient  of  the  best  least-squares  fit  to  the  data  with  orthogonal  polynomials 
error  measure  of  the  distance  between  the  polynomial  approximation  and  the  data 
the  difference  between  the  measured  ice  velocity  and  the  polynomial  representation 
at  local  time 

the  measured  ice  velocity  at  local  time 

drag  force  on  the  ice  due  to  fluid  shear 

downstream  gravitational  force  on  the  ice  due  to  its  weight 

acceleration  reaction  or  force  on  the  ice  due  to  the  apparent  or  added  mass 

force  on  the  ice  from  bank  resistance 

acceleration  due  to  gravity 

scalar  product  between  arbitrary  functions  g  and  h 

integrals  that  are  evaluated  numerically 

mass  of  ice  in  a  control  volume 

series  of  mutually  orthogonal  polynomials 

(th  polynomial  in  the  series  of  orthogonal  polynomials 

hydraulic  radius  associated  with  the  ice  passing  the  control  volume 

scalar  product  of  polynomial  P,  with  itself 

hydraulic  gradient  of  the  river 

global  time  for  the  ice  motion  in  a  river  reach 

local  time  at  a  point  measured  from  the  arrival  of  the  breaking  front  that  initiates 
ice  motion 

correction  to  7  at  each  step  of  the  iterative  solution 
dummy  time  variable  of  integration 
ice  thickness 
flow  velocity 

local  ice  velocity  at  a  fixed  location  represented  by  the  composite  polynomial 

global  ice  velocity  for  the  river  reach 

dummy  ice  velocity  variable  of  integration 

position  along  the  river  reach  in  a  fixed  coordinate  system 

initial,  final  position  of  the  ice  block 

location  of  the  breaking  front 

position  of  the  moving  coordinate  system  with  positive  distance  measured 

upstream  from  the  breaking  front 

distance  between  the  centers  of  neighboring  ice  blocks 

weight  assigned  to  the  measured  ice  velocity  at  local  time 

density  of  water 

density  of  ice 

bank  stress  on  the  ice  in  the  control  volume 
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Analysis  of  River  Ice  Motion  Near  A  Breaking  Front 

MICHAEL  G.  FERRICK,  PATRICIA  B.  WEYRICK.  SUSAN  T.  HUNNEWELL 


INTRODUCTION 

The  characteristics  of  river  ice  breakup  can  vary 
greatly  at  a  given  site  in  different  years,  and  between 
sites.  Damages  and  flooding  caused  by  breakup  generally 
increase  with  river  discharge  and  the  thickness  and 
competence  of  the  ice  cover.  As  the  force  of  the  river 
and  the  resistance  of  the  ice  cover  increase,  the  processes 
of  breakup  become  increasingly  dynamic.  In  such 
"dynamic  breakups”  both  the  ice  motion  and  the  river 
flow  exhibit  large  and  rapid  changes  with  time.  The 
importance  of  unsteady  flow  in  dynamic  breakup  has 
been  noted  by  several  authors  (e.g.,  Beltaos  and 
Krishnappan  1982,  Billfalk  1982,  Doyle  and  Andres 
1 979,  Ferrick  and  Mulherin  1 989,  Henderson  and  Gerard 
1981,Prowseetal.  1986,  Williamson  1989).  However, 
a  theory  of  dynamic  breakup  that  quantifies  the  time 
varying  flow  and  ice  conditions  is  not  available.  An 
essential  part  of  such  a  theory  would  describe  conditions 
near  the  front  of  the  breakup. 

The  following  description  and  analysis  of  dynamic 
ice  breakup  follows  Ferrick  and  Mulherin  (1989).  The 
“breaking  front”  in  a  dynamic  breakup  is  the  boundary 
between  intact  stationary  ice  and  moving  ice.  The 
breaking  front  travels  downstream,  and  when  it  passes 
a  given  point,  the  local  ice  cover  that  was  initially  at  rest 
also  begins  to  move  downstream.  The  breaking  front 
represents  a  discontinuity  in  the  ice  acceleration.  Ahead 
of  the  front  both  the  ice  velocity  and  acceleration  are 
zero,  and  behind  it  they  are  both  positive.  The  spieed  of 
the  breaking  front,  determined  by  the  force-resistance 
balance,  is  a  function  of  river  geometry  and  flow  and  ice 
characteristics.  Behind  this  front  the  ice  can  move  as  a 
single  sheet,  as  large  plates,  or  as  ice  rubble.  A  second 
front  called  the  “rubble  front”  separates  the  ice  sheet  or 
plates  from  the  brash  ice.  The  moving  ice  rubble  behind 
the  front  extends  above  the  downstream  ice  as  a  result 
of  its  thickness.  The  rubble  front  location  represents  a 
region  of  ice  convergence  and  is  readily  apparent. 

The  characteristics  of  dynamic  ice  breakup  depend 
on  the  mode  of  failure  of  the  ice  cover  A  failure  at  the 
supports  of  an  ice  cover  produces  a  sudden  bank-to- 


bank  release  of  the  ice.  This  ice  release,  which  we  call 
“support-dominated  dynamic  breakup,”  travels  rapidly 
downstream  at  a  speed  that  is  limited  by  the  dynamic 
wave  celerity,  but  is  much  greater  than  the  flow  veloc¬ 
ity.  Distinct  breaking  and  rubble  fronts  accompany  this 
type  of  failure,  with  the  rubble  front  following  the 
breaking  front  downstream.  The  distance  between  these 
fronts  varies,  depending  on  their  relative  instantaneous 
speeds.  The  breaking  front  in  a  support-dominated 
breakup  with  large  plate  or  sheet  movement  is  difficult 
to  observe,  a  consequence  of  minimal  ice  ridging  near 
the  front.  Conversely,  in  what  we  call  “strength- 
dominated  dynamic  breakup,”  the  brash  ice  from  breakup 
upstream  is  directly  involved  in  breaking  the  intact  and 
stationary  ice  cover.  The  ice  generally  fails  in  bending 
caused  by  forces  resulting  from  the  flow  and  from  the 
brash  ice  movement.  The  rubble  and  breaking  fronts  are 
coincident,  and  the  speed  of  this  combined  front  is 
usually  lower  than  the  flow  velocity.  However,  transi¬ 
tions  between  these  behaviors  are  common  during  a 
breakup,  allowing  plate  collisions  and  size  reduction 
ahead  of  the  ice  rubble  to  increase  the  apparent  speed  of 
the  nibble  front. 

The  motion  of  the  ice  near  the  breaking  front  during 
a  dynamic  ice  breakup  has  not  received  sufficient  at¬ 
tention.  This  motion  and  the  parameters  that  affect  it  are 
critical  for  development  or  assessment  of  theories  of 
dynamic  ice  breakup.  In  this  report  we  develop  an 
analysis  of  this  motion  and  demonstrate  the  results 
using  observed  data  from  a  controlled  dynamic  breakup 
of  the  Connecticut  River.  The  known  flow  conditions 
and  almost  rigid  body  motion  of  the  ice  sheet  near  the 
front  provide  a  relatively  simple  support-dominated 
breakup  for  study.  We  used  videotape  of  the  ice  motion 
to  identify  the  breaking  front  and  to  obtain  a  time  series 
of  ice  velocity  data  at  a  cross  section  of  the  river.  These 
data  typify  the  time  varying  ice  velocities  that  occur 
immediately  upstream  of  the  breaking  front.  An  empiri¬ 
cal  expression  for  ice  velocity  was  obtained  that  mini¬ 
mized  the  error  and  noise  in  the  data  and  identified  the 
primary  ice  motion.  This  methodology  should  be  gen¬ 
erally  applicable  to  any  dynamic  breakup. 


Our  analysis  of  the  ice  motion  relies  on  assumptions 
that  extend  these  data  at  a  point  to  a  complete  descrip¬ 
tion  of  the  motion  in  the  river  reach.  Many  of  the  results 
are  sensitive  to  the  speed  of  the  breaking  front,  a 
parameter  that  was  not  measured  locally.  However,  the 
speed  of  the  front  was  varied  over  a  wide  range  to 
explore  this  sensitivity.  The  time  variations  of  several 
parameters  related  to  the  ice  motion  were  obtained  for 
the  camera  site.  The  variability  of  the  ice  velocity  with 
time  indicated  a  need  to  formulate  and  study  an  equation 
of  motion  for  a  control  volume  at  this  site.  In  order  to 
more  completely  assess  conditions  near  the  front  we 
also  obtain  the  spatial  variations  of  several  parameters 
of  the  ice  motion  from  the  perspective  of  an  observer 
moving  with  the  breaking  front.  The  analysis  quantifies 
several  concepts  contained  in  the  largely  descriptive 
frameworkof  icebreakupof  Fenrickand  Mulherin  (1989). 

DATA  ACQUISITION 

A  controlled  ice  breakup  experiment  was  conducted 
on  the  Connect icu  t  Ri  ver  below  the  W ilder  Dam  (Ferrick 
andMulherin  1989)on  15-17March  1989.  A  controlled 
flow  surge  was  released  from  the  dam  that  caused  the  ice 
breakup  to  progress  downstream  with  the  front  of  the 
surge.  During  this  experiment  river  stage  and  ice  motion 
were  recorded  at  several  locations.  Video  cameras  were 
positioned  along  the  river  bank  to  record  the  ice  motion. 


Both  stage  and  ice  motion  data  were  collected  simulta¬ 
neously  at  a  site  located  26  km  downstream  of  the  dam. 
The  camera  at  the  26  km  observation  site  remained 
focused  on  a  position  in  the  river  for  almost  300  s 
starting  with  the  arrival  of  the  breaking  front.  This 
segment  of  the  tape  was  viewed  repeatedly  until  the 
paths  of  ice  motion  were  identified  and  drawn  on  a  clear 
sheet  overlying  the  video  monitor.  Points  on  the  ice  to 
be  used  for  velocity  measurement  were  chosen  by  their 
contrast  from  the  surrounding  ice.  Start  and  stop  points 
were  selected  on  each  ice  flow  line  near  the  center  of  the 
field  of  view  to  define  a  measurement  section.  The 
spacing  between  these  points  was  chosen  to  ensure  a 
sufficient  time  of  travel  for  precise  and  representative 
velocity  measurements,  and  a  length  scale  for  each  flow 
line  was  obtained  from  known  local  distances.  The 
videotape  was  window  dubbed  with  a  continuous  on¬ 
screen  digital  time  marker  to  provide  accurate  and 
consistent  travel  times. 

The  length  scale  divided  by  the  travel  time  through 
the  test  section  provided  the  estimated  ice  velocity  at  the 
mean  time  of  the  motion  (Fig.  1  ).The  data  obtained  thus 
represent  the  ice  velocity  as  a  function  of  time  at  one 
longitudinal  location  in  the  river.  Maximum  measure¬ 
ment  error  bounds  were  0.5  m  in  length  and  0.5  s  in 
elapsed  time.  Combining  these  errors  unfavorably  yields 
a  maximum  ice  velocity  error  of  slightly  less  than  7%  at 
the  peak  velocity,  corresponding  to  an  absolute  error  of 
0.08  m/s. 


Figure  F  Icc  velocity  data  from  the  controlled  breakup  of  the  Connecticut  River  in  March  1989  and  the  5th-order 
polynomial  description  of  the  data  obtained  from  a  least  squares  fit  of  the  sequence  of  orthogonal  polynomials. 
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LOCAL  ICE  VELOCITY  POLYNOMIAL 

We  will  now  seek  a  continuous,  differentiable  func¬ 
tion  to  represent  the  ice  velocity  \’i  data  that  were  ob¬ 
tained  at  the  measurement  section.  A  discrete  time  /„ 
measured  from  the  time  of  breaking  front  arrival  is 
associated  with  each  ice  velocity  measurement.  These 
data  are  represented  indicating  the  measured  ice 
velocity  of  th'e  nth  time  in  a  sequence  with  increasing 
time.  We  assume  that  the  measured  data  contain  a 
slowly  varying  component  that  provides  important  in¬ 
formation  about  the  changes  in  ice  velocity  with  time, 
and  a  comparatively  rapid  variation  of  relatively  small 
amplitude  that  represents  the  error  or  noise  in  the  data. 
The  task  is  to  approximateyffp)  with  a  function  V|(r)that 
contains  most  of  the  information  in  the  data  and  little  of 
the  error  or  noise.  In  our  analysis,  we  used  orthogonal 
polynomials  to  obtain  V|(r).  Details  of  the  application  of 
orthogonal  polynomials  to  fitting  the  ice  velocity  data 
can  be  found  in  the  appendix  to  this  report,  and  addi¬ 
tional  background  is  given  by  Conte  and  deBoor  ( 1 980). 

The  fitted  polynomial  for  ice  velocity  is  obtained  as 
a  linear  combination  of  orthogonal  polynomials.  We 
rewrite  this  expression  as  a  “composite”  of  the  orthogo¬ 
nal  polynomials  in  the  form 

>’iU  j  =  oo  +  +  ott  +  a^r  +.  .  . 

(1) 

where  T  is  local  time  measured  from  the  anrival  of  the 
breaking  front,  and  Aq  =  0  because  the  ice  velocity  is 
zero  at  7  =  0.  The  error  between  the  composite  polyno¬ 
mial  and  the  data  is  represented  by  a  time  series  for  each 
level  of  the  approximation.  The  error  between  the 
composite  polynomial  and  the  data  at  any  time  is 

ERROR  (0  =/(r„)- v,(tn)  (2) 


An  appropriate  degree  polynomial  can  be  found  by 
observing  both  the  absolute  magnitude  and  randomness 
of  the  error.  Polynomials  of  increasing  degree  incorpo¬ 
rate  additional  information  from  the  ice  velocity  data 
until  the  magnitude  of  the  maximum  error  ceases  to 
decrease  and  the  structure  of  the  error  is  random. 

The  ice  velocity  data  were  fitted  with  sets  of  orthogo¬ 
nal  polynomials  of  increasing  degree.  Based  on  the 
error  behavior,  the  set  through  degree  5  was  chosen  and 
a  plot  of  the  composite  polynomial  is  presented  with  the 
data  in  Figure  1 .  The  structure  of  the  error  is  shown  in 
Figure  2,  with  a  root-mean-square  (rms)  of  0.03  m/s. 
Errors  of  this  magnitude  are  in  the  expected  range  of 
measurement  error.  Higherdegree  polynomials  reduced 
the  rms  error  slightly,  but  did  not  reduce  the  maximum 
error  nor  significantly  increase  the  random  nature  of  the 
error. 


ICE  VELOCITY  AND  ACCELERATION 

During  dynamic  river  ice  breakup  the  ice  motion 
varies  in  both  time  and  space.  Data  are  not  available  to 
quantify  the  spatial  changes  in  ice  velocity  near  the 
measurement  location.  In  this  section  we  make  assump¬ 
tions  concerning  the  ice  motion  that  compensate  for  the 
Jack  of  multiple  ice  observation  points  in  the  reach  and 
provide  a  complete  description  of  the  motion.  We  relate 
the  local  time  and  ice  velocity  data  represented  by  the 
polynomial  to  a  “global”  reach  time  scale,  and  obtain  ice 
velocity  and  acceleration  for  the  reach  as  functions  of 
both  time  and  space. 

In  the  controlled  breakup  of  the  Connecticut  River, 
the  breaking  front  stalled  briefly  2  km  upstream  of  the 
observation  site  and  then  moved  quickly  through  the 
reach.  The  flow  disturbance  caused  by  the  ice  release 
upstream  initiated  the  breakup  of  the  study  reach.  We 


Figure  2 .  Structure  of  the  velocity  error  associated  with  the  5th-order  polynomial. 


will  make  the  following  assumptions  in  order  to  com¬ 
pletely  specify  the  ice  motion  near  the  breaking  front:  1 ) 
the  speed  of  the  breaking  front  is  analogous  to  water 
wave  celerity,  2)  this  breaking  front  celerity  was  con¬ 
stant  as  the  front  traveled  through  a  reach  of  about  2000- 
m  length  containing  the  observation  site,  and  3)  the 
measured  ice  velocity-time  behavior  was  the  same  at 
every  point  in  the  reach  when  offset  for  different  arrival 
times  of  the  breaking  front. 

With  these  assumptions  we  know  the  ice  velocity  as 
a  function  of  local  time  7  at  every  point  in  the  reach,  and 
local  times  are  related  by  the  relative  spatial  positions  of 
the  points  and  the  speed  of  the  breaking  front.  We 
arbitrarily  define  global  time  r  =  0  to  correspond  to  the 
breaking  front  at  location  j:  =  0.  If  Cj,  is  the  speed  of  the 
breaking  front,  then  the  front  arrival  at  any  location  x 
occurs  at  r  =  corresponding  to  7  =  0,  and  the  time 

scales  are  related  as 

7  =  /-;r/Cb  (3) 

In  differential  form  the  time  scales  given  in  eq  3  are 
related  as 


±=  1 

dl  Cb 


(4) 


where  V,  (x,t)  is  the  global  velocity. 

If  global  coordinates  x  and  t  are  specified,  we  ob¬ 


serve  that  the  corresponding  global  ice  velocity  Vi{jr,f) 
can  be  obtained  from  the  ice  velocity  polynomial  as 

V,(j:,r)  =  v,(7)  (5) 

Following  from  the  assumptions,  the  ice  velocity  is 
constant  along  parallel,  straight  lines  called  characteris¬ 
tics  in  the  x-t  plane  with  slope 


dx  _ 
dt 


Cb 


(6) 


shown  in  Figure  3.  Constant  ice  velocity  along  charac¬ 
teristics  given  by  eq  6  can  be  stated  as 


Cb-^1  V,  =  0 

dx 


(7) 


The  ice  velocity  polynomial  defines  the  motion  along  a 
vertical  line  in  Figure  3  that  extends  upward  in  time 
from  a  point  corresponding  to  the  observation  site.  A 
vertical  line  at  any  location  x  depicts  the  same  time  his¬ 
tory  of  ice  velocity  following  breaking  front  arrival.  The 
assumptions  have  expanded  the  description  of  motion  at 
a  point  contained  in  the  data  to  include  the  ice  motion  at 
all  points  in  an  x-t  plane  drawn  for  the  local  reach. 

Ahead  of  the  breaking  front  the  ice  is  stationary  with 
zero  acceleration,  while  behind  it  the  ice  undergoes 
positive  acceleration  and  begins  moving  downstream. 


Figure  3.  Sketch  of  the  x-t  plane 
showing  selected  lines  of  constant 
ice  velocity  including  the  breaking 
front.  Ice  velocities  at  fixed  loca¬ 
tions  are  the  same  through  time, 
offset  by  the  lag  in  breaking  front 
arrival.  Ice  velocities  with  distance 
from  this  front  are  the  same  at  any 
time. 
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The  ice  acceleration  Ai  is  obtained  by  differentiating 
eq  5 


A,  =  dT-  } 


cit 


dt 


dt 


vi  (0 

dt 


(8) 


which  can  be  evaluated  using  the  ice  velocity  polyno¬ 
mial.  The  derivative  of  the  ice  velocity  polynomial  with 
respect  to  7  is  the  local  ice  acceleration  as  a  function  of 
time  following  the  arrival  of  the  breaking  front.  Alter¬ 
natively,  by  expanding  dV^jdt  and  using  eq  7  we  obtain 


dt  dt  dx  3r  \  Cbl 


(9) 


With  eq  5  and  the  observation  that  the  local  acceleration 
dVi/dt  =  dvj(  7)/d7 ,  eq  8  and  9  are  equivalent. 

Although  the  speed  of  the  breaking  front  near  the 
observation  site  was  not  obtained  during  the  field  ex¬ 
periment,  the  average  speed  of  the  breaking  front  ob¬ 
served  over  a  1 2-km  reach  was  2.7  m/s,  which  included 
temporary  stalling  at  a  few  locations.  The  speed  of  the 
downstream  traveling  dynamic  wave  near  the  observa¬ 
tion  site  was  7.5  m/s.  The  breaking  front  speed  can  be 
nondimensionalized  by  the  flow  velocity  at  the  time  of 
breakup,  and  a  useful  range  of  dimensionless  front 
speeds  is  from  I ,  corresponding  to  the  flow  velocity,  to 
5,  corresponding  to  the  downstream  dynamic  wave 
speed.  The  total  ice  acceleration  at  every  point  in  the 
reach  is  presented  as  a  function  of  time  in  Figure  4  for 
dimensionless  breaking  front  speeds  of  2  and  5.  The 


zero  time  in  this  figure  corresponds  to  the  local  arrival 
of  the  breaking  front.  The  behavior  of  these  curves  is  the 
same  except  that  the  absolute  value  of  the  total  accelera¬ 
tion  increases  slightly  as  increases.  The  maximum 
acceleration  occurred  with  the  initial  motion,  deci  eas¬ 
ing  rapidly  in  the  first  minute.  For  the  remainder  of  the 
record,  acceleration  oscillated  about  zero. 


EQUATION  OF  MOTION  FOR  ICE 

We  will  now  write  a  simplified  equation  of  motion 
for  the  ice  cover  that  can  be  solved  analytically,  and 
begin  to  quantify  the  effects  of  various  parameters  on 
the  ice  motion.  In  this  initial  study  we  consider  the 
forces  on  the  ice  due  to  the  flow  but,  for  simplicity,  we 
neglect  the  response  of  the  flow  to  the  ice  motion.  The 
variations  in  flow  parameters  caused  by  ice  breakup  will 
be  quantified  in  future  work. 

The  flow  discharge,  depth  and  velocity  were  all 
increasing  through  the  period  of  ice  velocity  measure¬ 
ment.  The  rate  of  increase  of  river  stage  was  constant 
through  the  period,  but  the  increase  in  the  flow  depth 
was  only  a  few  percent  of  the  initial  value.  We  assume 
that  during  the  measurement  period  the  water  surface 
profile  of  the  surge  that  caused  the  breakup  moved 
downstream  as  a  monoclinal  wave  with  constant  speed 
and  unchanging  shape.  This  constant  rate  of  stage 
increase  thus  produces  a  constant  energy  gradient.  W ith 
constant  depth  and  energy  gradient  our  calculations 
indicate  an  increase  in  the  flow  velocity  of  about  25%  of 
the  initial  value  during  the  measurement  period.  Ihe 


F  igure  4.  Total  ice  acceleration  at  the  observation  site  for  assumed  dimensionless  bre  king  front  speeds  of  2  and  5. 


5 


increa'^e  in  flow  velocity  is  in  response  to  diminished 
resistance  :o  flow  as  the  ice  accelerates.  However,  the 
assumption  of  constant  velocity  at  the  mean  value  will 
be  made  to  siinplify  the  analysis  without  introducing  a 
large  error. 

During  tne  breakup  experiment  on  the  Connecticut 
River  the  sheet  broke  apart  downstream  of  the  measure¬ 
ment  location  and  provided  negligible  resistance  to  the 
moving  ice  upstream.  When  the  ice  run  finally  stalled, 
the  river  was  open  downstream  of  the  observation  site. 
Also,  force  dissipation  internal  to  the  ice  sheet  caused 
by  ridging  away  from  the  banks  was  not  common  near 
the  breaking  front.  Therefore,  we  can  neglect  down¬ 
stream  support  and  internal  force  dissipation  and  at- 
tribufe  all  ice  motion  resistance  to  the  river  banks. 
Because  of  daily  and  seasonal  changes  in  i.,er  stage, 
the  ice  sheet  exhibited  well-developed  shore  cracks 
near  each  bank  prior  to  the  controlled  brea  tup.  These 
cracks  were  roughly  parallel  to  the  banks,  but  highly 
irregular.  There  was  no  connection  between  the  cover 
and  the  banks  at  some  locations,  and  rigid  interlocking 
with  shorefast  ice  at  others.  To  obtain  a  tractable  equa¬ 
tion  of  motion  we  neglect  the  forces  transmitted  through 
the  cover  in  the  lorigitudinal  direction,  and  consider  the 
bank  resistance  term  as  typical  or  locally  averaged. 
Following  the  arrival  of  the  breaking  front,  the  ice 
fractured  into  small  pieces  near  the  banks  to  allow  the 
primary  sheet  or  floes  to  pass.  Tlie  amount  of  fracturing 
depended  on  the  size  of  these  large  ice  forms  relative  to 
the  local  river  width,  and  appeared  to  vary  with  ti  ne. 
These  processes  represented  an  energy  loss  to  the  mov¬ 
ing  ice  that  we  will  lump  into  the  bank  resistance. 

Consider  a  stationary  control  volume  that  spans  the 
river  cross  section  at  the  measurement  location  and  has 
a  unit  length  in  the  longitudinal  direction.  The  reason  for 
this  choice  of  width  is  that  the  sizes  of  the  ice  sheet  or 
plates  were  on  the  order  of  the  river  width.  The  unit 
length  is  arbitrary,  but  consistent  with  the  typical  or 
averaged  force  balance.  The  equation  of  motion  for  the 
ice  in  this  control  volume  can  be  written  as 

F’m  +  ^d+F'g-F'r='”l—  (10) 

dt 

where 

=  acceleration  reaction  or  force  due  to  the 
apparent  or  added  mass 
Fj  =  drag  force  on  the  ice  due  to  fluid  shear 

=  downstreamgravitationalfoicc.'uetotheweight 
of  the  ice 

F^  =  bank  resistance  force 

W|  -  mass  of  ice  in  the  control  volume. 

The  force  due  to  apparent  mass  can  be  parameterized  as 


where  is  an  apparent  mass  coefficient  to  be  deter¬ 
mined  experimentally  or  approximated  analytically. 
Wake  et  al.  ( 1 987)  performed  a  set  of  experiments  on  ice 
floes  with  various  aspect  ratios  and  found  =  0. 1  for 
floes  that  are  large  relative  to  their  thickness. 

'1  he  drag  force  on  the  ice  due  to  fluid  shear  will  be 
parameterized  using  a  hydraulic  radius.  An  idealized 
vertical  profile  of  water  velocity  in  an  ice  covered 
channel  exhibits  a  zero  velocity  at  each  boundary  and  a 
maximum  velocity  near  mid-depth.  The  precise  loca¬ 
tion  of  the  maximum  depends  on  the  relative  roughnesses 
of  the  cover  anJ  the  bed.  Physir''''y,  the  hydraulic 
radius  associated  with  the  ice  cover  lo  fhe  distance  from 
the  ice  to  the  velocity  maximum,  and  the  remainder  of 
the  depth  is  the  hydraulic  radius  associated  with  the  bed. 
The  hydraulic  radius  associated  with  the  ice  /?i(.v,.')  is 
obtaineo  from  the  Chezy  equation  as 

V'l(-riO)IV^-V^lUv)l  ('2) 
5f 


where 

C*|  =  dimensionless  Chezy  conveyance  coefficient 
associated  with  the  ice 
g  =  acceleration  due  to  gravity 
V  =  flow  velocity 
5f  =  hydraulic  gradient. 

The  drag  exerted  by  the  cover  on  the  flow  is  diminished 
during  breakup  when  >  0. 

With  the  parameterization  of  Ferrick  and  Mulherin 
( 1 989)  for  barik  resistance,  we  can  now  re  write  eq  ’  0  as 


pg  SSf 


/?,-(-  2lf, 
P  j 


ll  -  2  fiTb  =(l  -I- Cn,)  ^  (13) 

dr 


where 

p,pl  =  density  of  water  and  ice,  respectively 
B  =  river  width 
/|  =  ice  thickness 
m,  =  riBpi 

Tf,  =  bank  stress  on  the  ice. 


WTien  the  resistance  term  in  eq  13  is  small,  the  down¬ 
stream  component  of  the  weight  of  the  ice  can  cause  X'l  to 
exceed  V,  indicating  a  negative  /?|  and  downstream  drag 
of  the  ice  on  the  water  instead  of  resistance  to  the  flow. 

The  physical  parameters  p,  p,,  B.  C.|,  t,  and  the  hy¬ 
draulic  parameters  and  V  are  assumed  constant  at 
their  mean  values  for  the  initial  300  s  of  ice  motion.  The 
density  of  water  at  0°C  was  known,  and  river  width,  ice 
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Figure  5.  Hydraulic  radius  associated  with  the  ice  following  the  passage  of  the  breaking  front.  The  points  are  given 
by  R,  =  0.19  +  0.085  sin[n  +  (2n/I85)  (t  -  62)]. 


thickness  and  ice  density  were  measured  immediately 
prior  to  the  experiment.  These  measured  parameters 
displayed  minor  variability  relative  to  their  mean  values 
through  the  reach.  The  ice,  bed  and  combined  channel 
conveyances  for  the  reach  were  obtained  by  calibrating 
the  model  of  Ferrick  and  Mulherin  ( 1 989)  against  data 
for  both  open  water  and  ice  covered  flow  conditions. 
The  mean  energy  gradient  and  flow  velocity  values 
were  estimated  using  both  model  calculations  and  field 
measurements. 

The  hydraulic  radius  associated  with  the  ice  passing 
the  control  volume  was  found  from  eq  12  and  is 
presented  in  Figure  5.  f?i  diminished  dramatically  dur¬ 
ing  the  first  minute  of  ice  motion,  and  then  slowly 
oscillated  in  the  0. 1  -  to  0.3-m  range  for  the  remainder  of 
the  record.  The  rapid  decrease  in  f?|  for  the  control  vol¬ 
ume  demonstrates  that  the  flow  is  affected  by  the  ice 
motion,  and  as  mentioned,  further  work  that  quantifies 
this  coupling  at  breakup  is  needed.  The  resistance  to  ice 
motion  prevented  the  ice  velocity  from  reaching  the 
flow  velocity,  and  Fj  did  not  go  to  zero.  The  oscillating 
segment  of  the  data  is  closely  approximated  by  a  sine 
function  about  a  mean  of  0. 19  m,  with  an  amplitude  of 
0.085  m  and  a  period  of  1 85  s.  The  lack  of  any  attenu¬ 
ation  over  more  than  a  complete  period  indicates  either 


an  underdamped  response  of  the  ice  to  a  disturbance  or 
continuous  forcing  of  the  ice  velocity  oscillation.  A  free 
vibration  of  a  system  is  the  result  of  a  temporary  initial 
excitation  that  is  then  attenuated  by  damping  in  the 
system.  Alternatively,  the  harmonic  oscillations  in  a 
forced  vibration  are  maintained  by  the  forcing  function. 
We  do  not  expect  that  any  of  the  variables  taken  as 
constants  caused  the  periodic  behavior. 

TEMPORAL  SOLUTIONS  FOR  ICE  VELOCITY 
AND  BANK  STRESS 

In  this  section  we  will  address  the  question  of  whether 
the  oscillations  in  the  ice  velocity  indicate  a  free  or  a 
forced  vibration  by  developing  solutions  of  the  equation 
of  motion  subject  to  the  constant  parameter  approxima¬ 
tions  at  mean  values.  Initially,  we  will  remove  the  time 
dependence  from  the  equation  by  taking  dVfdt  =  0,  and 
obtain  equilibrium  ice  velocity  as  a  function  of  bank 
resistance.  Analytical  solutions  of  the  complete  equa¬ 
tion  are  then  developed  for  initial  conditions  of  ice  at 
rest  and  arbitrary  positive  and  negative  disturbances 
from  an  equilibrium  condition.  None  of  these  solutions 
exhibit  ice  velocity  oscillations  with  time.  Finally,  we 
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solve  the  equation  of  motion  algebraically  for  bank 
resistance  as  a  function  of  time,  given  ice  velocity  and 
acceleration  from  the  measured  data  and  eq  5  and  eq  8. 
Continuously  changing  bank  resistance  was  the  most 
likely  cause  of  the  ice  velocity  oscillations  contained  in 
the  data. 

Equilibrium  ice  velocity  as  a  function  of  bank  stress 
is  presented  in  Figure  6  for  the  hydraulic  conditions 
present  at  the  time  of  ice  breakup  at  the  observation  site. 
At  a  bank  stress  of  0.391  kPa  the  equilibrium  ice 
velocity  is  equal  to  the  flow  velocity.  The  maximum 
bank  stress  of  1 .87  kPa  corresponds  to  an  ice  velocity  of 
zero,  and  zero  resistance  corresponds  to  a  maximum  ice 
velocity  of  2.33  m/s.  The  constant  slope  of  much  of  this 
curve  follows  because  the  bank  stress  is  a  negative 
quadratic  function  of  the  ice  velocity  until  Vj  ap¬ 
proaches  V. 

In  most  cases  of  sheet  and  plate  motion  during 
breakup,  the  ice  velocity  is  lower  than  the  flow  velocity, 
and  eq  1 3  can  be  rewritten  as 


cVf  +  hV I  +  a  =  ( I  +  Cjff)  tn  | 

dt 


(14) 


where 


a  = 


pBV~ 

Cl\ 


+  Pi  2  tfCb 


b=-2  pBT/cIi 

c  =  pB/Cli  =--^ 

2V 

q  =  4ac  -b^<0 

We  are  primarily  interested  in  the  solution  for  ice 
velocity  at  a  point  because  these  results  can  be  directly 
compared  to  the  ice  velocity  polynomial.  Therefore, 
and  dV^ldt  in  eq  14  are  replaced  with  vj(  7)  and 
dvi/dl  using  eq  5  and  8,  and  eq  14  simplifies  to  an 
ordinary  differential  equation  in  local  time?  as 

cvf  -i-fcvi  +  a  =  (l +Cni)rwi[l -2^1^  (15) 

'  ^b/  dt 

Note  that  local  time  7 = 0  can  be  chosen  to  coincide  with 
the  arrival  of  the  breaking  front,  or  it  can  be  given  an 
arbitrary  definition,  and  i’|(  7)  is  a  general  local  ice  ve¬ 
locity.  With  constant  physical  and  hydraulic  parameters 
at  their  mean  values  we  can  separate  variables  in  eq  1 5 
and  write 

■  (1  +  C„)mi  j  '(Cb-v,)*, 

— p.  -  -b — ('<» 

CV’I  +  -I- a 

where  t  and  (q  are  variables  of  integration,  and  the  limits 
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of  the  integrals  indicate  corresponding  local  times  and 
ice  velocities  at  the  control  volume. 

We  initially  integrate  eq  16  for  the  case  of  vj «  C5 
and  obtain 


( 1  +  Cm) 


wi 


_2 


tanh 


-1  2cv|  +  h 


v'l 


(17) 


where  vi^  is  the  ice  velocity  at  /  =  0.  After  some  ma¬ 
nipulation,  eq  17  can  be  inverted  to  obtain  V|  as  a 
function  of  time: 


I’l  = 

where 


^’lo  + 


-^1  ^ 
2c 


l-Ml) 


\  +AiA2 


(18) 


A,  =tanh  I _ Gj - 

\(1  +Cm)w,(-2) 


2c-V|o  +  h 


Integration  of  eq  16  without  simplification  yields 


This  equation  for  7  as  a  function  of  vj  is  not  readily  in¬ 
verted.  However,  the  solution  expressed  by  eq  19  can 
be  obtained  for  all  ice  velocities  except  the  equilibrium 
value.  At  equilibrium  the  ice  velocity  does  not  vary  with 
time,  and  the  one-to-one  correspondence  between  vj 
and?  no  longer  applies.  The  solution  of  the  equation  of 
motion  can  be  obtained  similarly  for  the  case  of  V'l  >  T. 

The  solutions  expressed  in  eq  18  and  eq  19  are 
presented  in  Figure  7  for  an  initial  ice  velocity  of  zero, 
several  values  of  bank  stress,  and  dimensionless  break¬ 
ing  front  speeds  of  2,  5  and  “.  As  the  dimensionless 
breaking  front  speed  decreases  the  approach  to  equilib¬ 
rium  occurs  more  rapidly  if  all  other  parameters  are  the 
same.  The  analytical  solution  closely  approximates  the 
ice  velocity  data  with  T|,  =  0.55  kPa,  except  for  the 


Figure  7.  Analytical  solution  for  the  increase  in  ice  velocity  from  rest  as  a  function  of  time  for  several 
values  of  bank  resistance  and  several  dimensionless  breaking  front  speeds. 


9 


Figure  8.  Analytical  solution  for  the  return  to  equilibrium  after  positive  and  negative  initial 
disturbances  for  a  single  bank  resistance  and  several  dimensionless  breaking  front  speeds. 


velocity  oscillations.  These  same  solutions  are  given  in 
Figure  8  for  this  value  of  bank  stress  and  four  arbitrary 
initial  ice  velocities  that  represent  positive  and  negative 
disturbances  from  equilibrium.  The  relative  behavior  of 
these  solutions  follows  that  of  the  curves  in  Figure  7. 

In  order  to  understand  these  results  we  consider  two 
breakup  conditions  with  the  same  local  velocity  and 
local  acceleration  behavior,  and  with  all  parameters 
equal  except  for  the  breaking  front  speed.  As  the  speed 
of  the  breaking  front  decreases,  approaching  the  ice 
velocity,  spatial  changes  in  the  ice  velocity  increase, 
and  the  total  acceleration  and  the  effective  inertia  of  the 
ice  are  reduced  (eq  1 5).  Under  these  conditions  equilib¬ 
rium  is  quickly  reestablished  following  a  disturbance, 
and  therefore,  this  equilibrium  is  very  stable.  Con¬ 
versely,  increasing  Cj,  corresponds  to  increasing  the 
inertia  of  the  ice,  and  increased  time  to  restore  equilib¬ 
rium.  There  are  no  oscillations  of  the  ice  velocity  with 
time  in  any  of  these  solutions. 

Numerical  model  results  of  Ferrick  and  Mulherin 
(1989)  indicate  that  the  bank  stress  increases  with  time 
in  response  to  the  arrival  of  a  flow  surge,  reaching  a 
maximum  at  the  arrival  of  the  breaking  front.  We  can 
develop  curves  of  bank  stress  with  time  beginning  at 
7=  0  for  the  observation  site  from  eq  15  ,  given  the  ice 
velocity  and  acceleration  relations  as  input.  The  results 
are  presented  in  Figure  9  for  apparent  mass  coefficient 


values  of  0.0  and  0. 1  and  dimensionless  breaking  front 
speeds  of  2  and  5,  covering  the  range  of  probable 
parameter  values.  These  curves  are  very  similar,  with 
the  greatest  differences  early  in  the  motion  when  the  ice 
acceleration  was  near  its  maximum  value.  The  mini¬ 
mum  bank  resistance  occurred  at  about  30  s  into  the 
motion.  After  this  initial  period  the  bank  shear  is  out  of 
phase  with  the  ice  acceleration  and  precedes  changes  in 
ice  velocity.  The  maximum  bank  stress  in  this  record 
was  about  38%  of  that  preceding  the  failure,  causing  ice 
deceleration  and  increasing /f|.  Continuously  changing 
bank  shear  due  to  variable  ice  plate  size  was  a  likely 
cause  of  variable  ice  velocity,  and  it  is  the  only  process 
in  this  model  that  can  account  for  this  variability. 

Due  to  the  river  control  operations,  the  discharge  and 
the  downstream  force  on  the  ice  gradually  diminished 
following  the  period  of  ice  velocity  measurement.  At  1 7 
minutes  after  the  arrival  of  the  breaking  front,  the  ice 
plate  motion  arrested  at  the  measurement  location.  The 
stall  occurred  at  the  peak  river  stage,  with  known  /?|  and 
5f  closely  approximated  by  the  bed  slope.  The  bank 
stress  at  the  stall  can  be  estimated  using  eq  1 3  as  1 .53 
kPa,  about  82%  of  that  prior  to  the  initial  motion.  A  large 
flow  release  on  the  following  day  yielded  a  slightly 
higher  stage  at  the  downstream  force  peak  with  T|,=  1 .59 
kPa,  and  the  large  plates  did  not  fail.  We  conclude  that 
large  hydraulic  forces  were  the  cause  of  the  ice  motion 
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Time  (s) 

Figure  9.  Bank  stress  following  the  passage  of  the  breaking  front  at  dimensionless  speeds  of  2  and  5.  both 
considering  and  neglecting  the  force  due  to  apparent  mass. 


and  not  simply  high  river  stage.  The  eventual  stall  of  the 
run  was  due  to  both  diminishing  forces  and  high  bank 
resistance. 


the  total  time  of  motion  t  must  be  related  to  7,  the  time 
since  the  passage  of  the  breaking  front  at  the  present 
location  of  the  block.  These  time  scales  are  related  in 
integral  form  from  eq  4  as 


SPATIAL  VARIATIONS  NEAR  THE 
BREAKING  FRONT 


dt 

Cb  -  Vi  (r) 


=  /i 


(20) 


We  will  now  examine  the  spatial  variation  of  several 
parameters  that  describe  the  ice  motion  near  the  break¬ 
ing  front.  An  analysis  of  the  spatial  variability  comple¬ 
ments  the  earlier  one  that  provided  temporal  changes, 
and  completes  the  description  of  the  ice  motion.  Spatial 
variations  of  motion  parameters  behind  the  breaking 
front  are  related  to  those  in  time  at  a  point.  In  Figure  3 
we  observe  that  the  ice  position,  velocity,  and  accelera¬ 
tion  upstream  of  the  breaking  front  become  constant  in 
time  when  spatial  position  is  relative  to  the  front.  The 
lack  of  ice  ridging  and  rubbling  near  the  front  in  the  field 
indicates  that  the  ice  sheet  was  not  undergoing  signifi¬ 
cant  ice  convergence.  We  find  that  ice  convergence 
local  to  the  breaking  front  depends  on  the  speed  of  the 
front,  and  note  the  implications  of  this  result  on  breakup 
behavior. 

In  order  to  obtain  the  state  of  motion  of  an  ice  block. 


where  r  =  0  is  defined  as  the  time  of  initial  motion  of  the 
block  and  t  is  a  dummy  variable  of  integration.  We 
observe  that  /  in  eq  20  is  known,  the  unknown  Tappears 
as  the  upper  limit  of  the  integral  and7<  r.  The  solution 
for?  can  be  obtained  by  trial,  starting  from  an  initial 
estimate?  =  (1  -  vi{t)IC^)t,  and  repeatedly  applying 
corrections  until  /i  =  f.  The  correction  to?  at  each  itera¬ 
tion  is 


At  = 


(/]  -t) 


(21) 


and  the  updated  value  is  obtained  as  ?  =  ?  -  A  ?.  The 
degree  of  the  polynomial  represented  by  V|(/ )  in  eq  20 
will  depend  on  the  data.  Therefore,  we  evaluate  /|  using 
Gaussian  quadrature  with  limits  changed  to  the  standard 
interval  (-1,1): 
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function  of  the  time  of  block  motion  t  fora  range  of  dimensionless  breakingfront 
speeds. 


The  timeTas  a  function  of  the  time  of  block  motion  t  is 
given  in  Figure  1 0  for  the  Connecticut  River  ice  velocity 
polynomial  and  a  range  of  dimensionless  breaking  front 
speeds.  At  small  values  of  time  and  for  large  values  of 
CyJV,  the  two  time  scales  are  nearly  the  same.  However, 
for  small  dimensionless  front  speeds,  these  time  scales 
diverge  with  increasing  time. 

The  position  of  an  ice  block  that  has  been  in  motion 
for  time  t  can  be  obtained  by  relating  position  and 
velocity  and  then  changing  the  variable  of  integration 
with  eq  4  and  5 

=  I'  V,  (xj)  dt  =  Cb  r  di  =  h  (23) 

Cb  -  v|  (t) 


Given  t  the  corresponding  value  of  7  is  known  from  eq 
20-22,  and  the  limits  of  integration  in  h  are  specified. 
The  final  block  position  .Vf  can  then  be  obtained  explic¬ 
itly  given  the  initial  position  x^.  As  before,  I2  is  evalu¬ 
ated  using  Gaussian  quadrature  on  the  standard  interval 
(-1,1).  The  final  position  of  a  block  relative  to  the 
breaking  front  can  be  determined  from  the  time  of 
motion  and  speed  of  the  front  together  with  the  final 
position  of  the  block,  or  from  the  speed  of  the  front  and 
r  at  the  final  position  of  the  block.  The  velocity  of  this 
ice  block  follows  from  eq  1 ,  and  block  acceleration  is 
obtained  from  eq  8  . 

The  river  ice  sheet  is  idealized  as  a  collection  of 
independent  blocks  initially  at  rest  with  known  fixed 
coordinate  positions  x,  and  having  a  uniform  initial 
spacing  &x^.  A  moving  coordinate  system  traveling 
with  the  breaking  front  passes  x  =  0  at  global  time  t  =  0. 
The  moving  distance  coordinate  is  positive  upstream 
from  the  breaking  front,  and  the  fixed  and  moving 
systems  are  related  as 
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'bf  -'m  (24) 

where  the  breaking  from  location  in  the  fixed  system 
=  Cfjf.  Comparing  eq  24  with  eq  3  we  observe  that 
=  C|,7,  where  7  is  the  local  time  at  position  x.  Given  a 
dimensionless  breaking  front  speed  and  the  local  time  7 
at  the  final  position  of  a  block,  the  corresponding 
distance  from  the  front  is  known  and  a  spatial  view  of 
any  parameter  can  be  obtained.  We  will  track  many 
blocks  from  stationary  initial  positions  to  a  location  and 
a  state  of  motion  at  a  later  time.  The  results  obtained 
below  are  all  presented  with  distance  .r^,  from  the 
breaking  front. 

The  time  of  block  motion  as  a  function  of  distance 
from  the  breaking  front  is  presented  in  Figure  11  for  a 
range  of  dimensionless  speeds  of  the  front.  As  the  front 


speed  approaches  the  flow  velocity,  the  time  of  motion 
rapidly  increases  immediately  behind  the  breaking  front. 
Blocks  that  are  in  motion  for  a  long  time  experience 
many  collisions  and  typical  sizes  are  greatly  reduced. 
As  the  speed  of  the  front  increases,  the  time  of  block 
motion  at  a  comparable  distance  from  the  front  dimin¬ 
ishes,  allowing  fewer  opportunities  for  collisions  and 
resulting  in  larger  typical  sizes.  A  block  that  has  been  in 
motion  for  100  s  behind  a  front  with  a  dimensionless 
speed  of  5  is  about  8  times  farther  from  the  front  than  the 
same  block  in  a  breakup  with  a  front  speed  of  1 . 

Ice  velocity  and  acceleration  as  a  function  of  dis¬ 
tance  behind  the  breaking  front  are  given  in  Figures  1 2 
and  1 3,  respectively,  for  a  range  of  dimensionless  front 
speeds.  As  can  be  seen  in  Figure  12,  the  spatial  ice 
velocity  gradient  decreases  as  the  speed  of  the  breaking 


Figure  I  /  Time  of  block  motion  as  a  function  of  distance  behind  the  breaking  front  for  a  range  of 
dimensionless  speeds  of  the  front. 
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Distance  Upstream  from  Breaking  Front  (m) 


Figure  13.  Total  acceleration  as  a  function  of  distance  behind  the  breaking  front  for  a  range  of 
dimensionless  speeds  of  the  front. 


Figure  14.  Ice  convergence  as  a  function  of  distance  behind  the  breaking  front  for  a  range  of 
dimensionless  speeds  of  the  front. 


front  increases.  Paralleling  the  results  at  a  point,  the 
spatial  acceleration  relations  depicted  in  Figure  13 
indicate  that  the  maximum  acceleration  occurs  at  the 
breaking  front  for  all  front  speeds,  and  that  the  total 
acceleration  decreases  more  rapidly  behind  the  front  as 
(?(,  decreases. 

Ice  convergence  is  an  important  characteristic  that 
distinguishes  the  types  of  dynamic  breakup.  When  the 
breaking  front  passes  block  i  it  begins  to  accelerate.  The 
spacings  between  this  block  and  its  neighbors  AXj  _  ]  and 
Axj  then  begin  to  change.  We  will  define  a  quantitative 
measure  of  ice  convergence  as  the  final  block  spacing 
divided  by  the  initial  spacing  Ax/Axg.  During  breakup 
this  ratio  is  generally  <  1,  indicating  a  reduced  ice 
surface  area  relative  to  the  initial  sheet.  This  ice  con¬ 
vergence  measure  quantifies  the  extent  of  the  surface 
area  loss  that  is  necessary  for  ice  continuity. 

Ice  convergence,  like  the  other  parameters  of  ice 
motion,  also  develops  a  steady  state  in  a  coordinate 
system  moving  with  the  breaking  front.  As  before,  the 
results  presented  in  Figure  1 4  are  for  a  range  of  dimen¬ 
sionless  breaking  front  speeds.  Ice  convergence  occurs 
everywhere  in  the  moving  pack,  and  generally  increases 
with  distance  behind  the  front.  The  spatial  gradient  of 
ice  convergence  quantifies  the  local  loss  of  ice  surface 
area  per  unit  length.  The  maximum  spatial  gradient 
occurs  immediately  behind  the  breaking  front  in  all 


cases.  It  decreases  significantly  as  the  speed  of  the  front 
increases.  Only  when  C^fV  <  2  is  the  gradient  near  the 
front  much  greater  than  at  other  points  upstream,  and 
these  values  are  more  than  an  order  of  magnitude  larger 
than  those  at  C^/V^4.  At  dimensionless  breaking  front 
speeds  of  less  than  1 ,  Ax/Axp  can  be  negative  with  these 
ice  velocity  data,  indicating  that  the  relative  positions  of 
the  ice  blocks  are  interchanged.  The  combination  of 
large  travel  time  and  significant  ice  convergence  occurs 
at  low  dimensionless  breaking  front  speeds. 

On  the  Connecticut  River  the  ice  convergence  near 
the  front  was  small,  as  evidenced  by  minimal  ridging 
within  the  sheet.  These  observations  match  the  com¬ 
puted  results  if  the  local  breaking  front  speed  ap¬ 
proached  the  dynamic  wave  speed.  At  the  other  ex¬ 
treme,  as  the  front  speed  approaches  the  flow  velocity, 
the  ice  convergence  becomes  large  at  the  breaking  front 
as  is  typical  in  strength-dominated  breakup.  The  di¬ 
mensionless  speed  of  the  rubble  front  on  the  Con¬ 
necticut  River  averaged  just  under  1.5.  Transitions 
between  strength-  and  support-dominated  breakup 
provide  an  opportunity  for  plate  collisions  and  size 
reduction  to  occur  ahead  of  an  existing  rubble  front. 
Therefore,  the  rubble  front  speed  can  exceed  the  flow 
velocity  when  the  ice  velocity  does  not,  and  the  sepa¬ 
ration  distance  between  the  breaking  and  rubble  fronts 
is  limited. 
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CONCLUSIONS 

In  this  report  we  quantify  parameters  of  the  ice 
motion  near  the  breaking  front  for  dynamic  breakup  of 
a  river  ice  cover.  Dynamic  breakup  of  river  ice  is 
generally  accompanied  by  large  and  rapid  changes  in 
both  the  river  flow  and  the  ice  motion.  TTie  concept  of 
the  breaking  front  as  the  boundary  between  moving  and 
stationary  ice,  and  the  motion  of  the  ice  behind  it  have 
been  investigated  using  data  from  the  Connecticut  River. 
Our  analysis  required  the  assumptions  of  constant 
breaking  front  speed  through  a  short  river  reach  and 
identical  ice  velocity  behavior  at  each  point  in  the  reach. 

The  ice  motion  was  initially  studied  from  the  per¬ 
spective  of  a  stationary  observer  on  the  river  bank.  The 
ice  acceleration  was  discontinuous  with  the  passage  of 
the  breaking  front,  going  from  zero  to  the  peak  value, 
and  total  acceleration  was  related  to  local  acceleration 
by  the  parameter  (1  -  vi/Cb).  The  hydraulic  radius  as¬ 
sociated  with  the  ice  diminished  dramatically  during  the 
initial  60  s  of  the  ice  motion,  indicating  a  rapid  decrease 
in  the  flow  resistance  of  the  ice.  An  equation  of  motion 
for  the  ice  and  solutions  of  this  equation  with  constant 
parameters  were  developed.  Equilibrium  ice  velocity 
occurs  during  steady  motion  and  is  a  function  of  bank 
stress.  When  the  ice  motion  is  unsteady,  oscillations  in 
the  ice  velocity  do  not  occur.  The  observed  ice  velocity 
oscillations  require  forcing,  and  the  most  likely  cause 
of  this  forcing  is  variable  bank  resistance.  If  all  other 
parameters  are  the  same,  the  rate  of  approach  to  equilib¬ 
rium  increases  as  the  speed  of  the  breaking  front  de¬ 
creases. 

The  ice  motion  was  also  observed  from  the  breaking 
front  over  a  range  of  dimensionless  front  speeds  to 
quantify  the  spatial  variations  of  conditions  near  the 
front.  Immediately  behind  the  front  the  time  of  ice 
motion,  velocity,  acceleration  and  convergence  do  not 
vary  with  time  but  are  strongly  related  to  the  speed  of  the 
front.  As  the  speed  of  the  breaking  front  decreases,  the 
time  of  ice  motion  at  a  given  distance  upstream  of  the 
front  increases.  Ice  convergence  occurred  everywhere 
in  the  moving  pack,  but  it  was  always  at  a  minimum 
immediately  adjacent  to  the  front.  As  the  speed  of  the 
front  approached  the  downstream  dynamic  wave  speed 
Cd,  the  ice  convergence  was  small  over  a  significant 
distance  upstream  of  the  front.  Minimal  ice  conver¬ 
gence  near  the  front  indicates  that  ice  continuity  can  be 
satisfied  without  the  development  of  ridges  or  other 
features  that  are  obvious  to  an  observer.  If  the  speed  of 


the  breaking  front  is  less  than  or  equal  to  the  ice  velocity, 
the  ice  convergence  near  the  front  is  large  and  the 
location  of  the  front  is  distinct. 

Minimal  ice  convergence,  short  travel  times  and 
large  block  sizes  near  the  breaking  front  are  typical  of 
support-dominated  breakup.  These  characteristics,  ob¬ 
served  at  the  ice  velocity  measurement  site,  are  consis¬ 
tent  with  this  analysis  for  Cb  =  Cd  »  Vi.  Conversely, 
this  analysis  indicates  that  if  Cb  =  V\  «  Cd,  the  ice 
convergence  and  time  of  ice  motion  are  large  and  block 
sizes  are  small  near  the  front,  characteristics  that  are 
typical  of  strength-dominated  breakup. 
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APPENDIX:  ORTHOGONAL  POLYNOMIALS 

In  the  following  discussion  we  summarize  a  more 
detailed  treatment  of  orthogonal  polynomials  by  Conte 
and  deBoor  ( 1 980)  and  apply  it  to  fitting  our  ice  velocity 
data.  Orthogonal  polynomials  offer  several  advantages 
for  fitting  data.  The  property  of  orthogonality  avoids  the 
condition  problem  in  the  normal  equations  that  yields 
poor  least-squares  approximations.  These  polynomials 
satisfy  a  recurrence  relation  between  the  next  higher 
degree  polynomial  and  the  two  previous  polynomials  in 
the  sequence,  allowing  progressive  development  with 
great  effic  iency .  The  error  between  the  composite  poly¬ 
nomial  and  the  data  is  represented  by  a  time  series  for 
each  level  of  the  approximation.  By  comparing  the 
amplitude  and  structure  of  the  error  it  is  relatively  easy 
to  assess  the  advantage  of  higher  degree  approxima¬ 
tions.  Finally,  these  polynomials  are  readily  integrated 
and  differentiated. 

The  ice  velocity  data  were  obtained  at  discrete  times 
in  the  interval  (0,  ^  positive  weighting  function 

n  fr)  is  defined  on  this  interval,  and  a  scalar  product 
between  the  arbitrary  functions  g  and  /i  can  be  defined 
as 

.V 

<  H-h>=  X  ^  h{ln)w{tn)  j 

where  r | ,  fs . 'n  times  on  the  interval.  The 

functions  and  are  orthogonal  on  this  in¬ 
terval  i  f  <g,/i>  =  0.  Given  that  the  polynomials  P  \  {t), 

Piit). ...  are  all  orthogonal  to  each  other  and  that  each 
P,  is  a  polynomial  of  exact  degree  i,  then  any  polynomial 
V|(/)  of  degree  <  k  can  be  written  as 

v,(f)  =  doP^it)  +  diP,(r)  +  . . .  +  (A2) 

with  coefficients  dg,  d, . uniquely  determined  by 

v,(r). 

In  our  analysis  we  will  obtain  a  composite  polyno¬ 
mial  v^{t)  that  is  equivalent  to  eq  A2  and  provides  the 
best  least  squares  approximation  of  the  ice  velocity 
data.  We  can  express  the  squared  error  in  the  polynomial 
approximation  as 

N 

<A‘)~  v^t)>  =  ^  l/(fn)- Vl(tn)fw(tn) 

n  =  I 

(A3) 

and  seek  a  polynomial  of  degree  k  that  minimizes  this 
error.  The  ice  velocity  at  r  =  0  is  well  known  [/(O)  =  0] 
compared  to  the  other  data.  Therefore,  we  specify  that 
all  the  data  receive  equal  weighting  and  require  that 


i’[(0)  be  sufficiently  close  to  zero  by  giving  that  point  a 
greater  weight, 

w(0) »  1 

~  ^  ®  —  ^max 

Larger  values  of  m'(0)  decrease  the  error  at  r  =  0,  but  the 
other  errors  generally  increase.  Substituting  eq  A2  into 
eq  A3  yields  the  error  £  as  a  function  of  the  coefficients 
dg, . . . .  dj.: 


£  (dg.  .  .  .  ,  dk)  =  <At)  -  dgPg  -  .  .  .  - 

~  ^0^ 0  ~  ~  k  ^  (A4) 


To  minimize  £  we  take  V£  =  0  and  obtain 
N 

Z  [/(fn)-vKfn)M^n)/’i=0 
n  =  1 

or  equivalently  because  the  polynomials  are  ortho¬ 
gonal, 

<£,,/>  =  <£;,v,>  =  </>i,£i>d,*  (A5) 

where  the  *  notation  indicates  the  best  coefficients. 
Rearranging  eq  A5  we  obtain 


</-!./> 

<  Pf,  P'l  >  Sj 


(A6) 


The  set  of  orthogonal  polynomials  is  developed  from 
the  three-term  recurrence  relation 


(A7) 


with 


and 


B,  =  ^  for  1  =  0,  1, 2 .  Jt  -  1 

Si 


Ci  =  ^i-for(  =  1,2 . k-\ 

Si  - 1 


Equations  A6  and  A7  require  that 


Si=  Z  [/'i  (fn)f  Mfn)  ^  0 
n  -  I 

To  initiate  the  polynomial  sequence  we  define 


(A8) 


P_iW  =  0 


Po(t)  =  1 

without  the  loss  of  generality. 
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